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(54) Method for determining the attenuation of a PCM signal over a digital channel 



(57) A method of determining digital channel atten- 
uation; comprising the steps of: receiving a known train- 
ing sequence of PCM codes, which PCM codes are 
subjected to the attenuation within the digital channel; 
quantizing the received known training sequence of 
PCM codes according to a predetermined thresholding 
procedure; identifying identical PCM codes created as a 
result of the thresholding procedure; and, determining 
the attenuation of the digital channel based upon the 
identification of identical PCM codes. A method is also 
disclosed for determining a digital channel PCM code 
transformation comprising receiving a known training 
sequence of PCM codes, which PCM codes are sub- 
jected to the PCM code transformation within the digital 
channel, quantizing the received known training 
sequence of PCM codes according to a predetermined 
thresholding procedure, and determining the transfor- 
mation of transmitted codes to those received . A 
method is also disclosed for improved echo cancellation 
in a communications network having an analog and a 
digital modem, comprising saving codes transmitted 
from the digital modem to the analog modem for echo 
cancellation, transforming, by a mapping table, codes 
transmitted from said digital modem to codes received 
by the analog modem, and, using the received codes as 
a reference signal for cancellation of echo. A method of 
improved spectral shaping using a transmit shaping 
transfer function in a communications network having 
an analog and a digital modem, comprising, transform- 
ing, by a mapping table, codes transmitted from the dig- 
ital modem to codes received by the analog modem, 



using the received codes for transformation to their lin- 
ear value equivalent representations, and, applying the 
linear value representations to the transmit shaping 
transfer function. 
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Description 

Field of the Invention 

s This invention pertains generally to modem technology, more specifically to PCM modem technology, and more 

specifically, to a method for the discovery of the digital attenuation from a set of received PCM samples corresponding 
to a known set of transmit PCM samples. 

Background of the Invention 

w 

This invention provides a method for the discovery of the digital attenuation from set of received pulse code modu- 
lation (PCM) samples corresponding to a known set of transmit PCM samples. TTie principal feature set used is the 
knowledge of which received PCM codes have become indistinguishable (i.e., are identical) as a result of the attenua- 
tion mapping. An alternative but equivalent feature set would be use of the absent received PCM codes. 
is The question of digital attenuation PAD mapping has been raised at recent meetings of the' Telephone Industry of 
America (TIA) PCM Modem Ad-hoc Committee. Many participants have been hoping for a single industry standard 
mapping solution. We provide a solution, which will allow PCM modems to function properly throughout the entire tele- 
phone network including with private telephone equipment (PBX's and key systems). 

20 TIA 464A 

TIA 464A is the industry standard plan for private telephone equipment. Among other recommendations, this spec- 
ification provides the recommended losses between T1 lines and other devices. If a PBX is designed according to the 
loss-plan specification as outlined in sections 4.8.4 and 4.8.5 of 464A, then it should not be necessary to attenuate the 
25 PCM on incoming T1 lines for on-premise (ONS) connections, since the recommended 3-dB loss is typically imple- 
mented in the analog circuitry of the ONS codec. 

However, Inter-Tel's experience with customers using T1 lines suggests that the recommended 3-dB of insertion 
loss is inadequate and that T1 attenuation pads are indeed necessary. While the standard -3 and -6 dB pads are usually 
adequate, Inter-Tel has provided a wider-range of insertion loss as options for PBX customers. Inter-Tel has imple- 
30 mented digital gain control on its T1 line cards using two different methods: 

ROM-based mapping 
DSP-based algorithm 

35 ROM Based Mapping 

Inter-Tel first implemented attention pads via an EPROM circuit, which used the incoming PCM code as 8 bits of an 
address and additional address bits selecting the attenuation in 1 dB increments. The output of the EPROM was the 
attenuated digital PCM code. In this implementation, the attenuation range covered from 0 dB to -12 dB in 1 dB steps. 
40 Experience with customers shows that 0 to -6 dB seems to be the most-often used range of T1 pad values. Taking 
into account the fixed 3-dB of loss included with our ONS card, the net insertion loss for to Digital CO connection is then 
-3dB to -9dB, which mimics well the typical net losses experienced by customers on ONS to Analog CO connections. 

All attenuation tables in the EPROM preserved the LSB of the PCM to preserve the state of the Robbed Bit Sign- 
aling (RBS), if present. However, the preservation of the LSB is unnecessary if no RBS information passes through the 
45 ROM-based gain control circuitry. 

The ROM-based look-up tables were generated by a custom C program using the following algorithm: 

1 . Expand incoming 8-bit,u-Law PCM to corresponding 14-bit, signed integer linear value, x 

2. Calculate output, y, based on equation: y=INT [x*10 (gain/20) )+0.500] 

so 3. Compress y back into its corresponding 8-bit u-Law value using G.71 1 decision values. 

In addition to performing u-Law to u-Law digital gain control (DGC) through look-up tables, the EPROM also per- 
forms DGC through look-up tables for the following compression schemes: 

55 • u-Law-to-A-law 
A- law- to -u-Law 
A-law-to-A-law 
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DSP Based Mapping 

With the advent of DSP's, for their greater flexibility, the latest implementation uses the Analog Devices 21xx fixed 
point family of processors. The PCM codes are converted to/from linear values by the ADSP-21xx internal companing 

5 hardware. Inter-Tel 's algorithm normalizes the linear value before applying the attenuation multiplier. The attenuation is 
selectable as a linear 1.15 fractional multiplier. The user interface selects the attenuation in fixed, 1 -dB steps over a sim- 
ilar range as the previous implementation. The user configured dB attenuation is converted to the linear fractional mul- 
tiplier by an algorithm. Moreover, the instructions in the Analog Devices DSP can use unbiased rounding, truncation 
and in newer family members, biased rounding. Inter-Tel's implementation uses the unbiased rounding option (RND) 

10 during the MAC instruction. This flexibility in rounding and truncation obviously can have a substantial effect on the 
actual attenuation mapping, with the additional caveat of having potentially different mappings for positive and negative 
PCM codes. 

Detailed Description of the Preferred Embodiment 

15 

Detection pf the Feature Set 

After a receiver's equalizer has been trained, a known transmit sequence of PCM codes can be sent and received 
using either ail or a useful subset of PCM codes. (This sequence may need to be sent multiple times to cover all 6 or 
20 12 bit positions in a robbed bit signaling system.) The receiver can save the linear values for each received PCM code. 
From the known step size between adjacent PCM codes, the receiver can determine those PCM samples that it 
received which are indistinguishable. The indistinguishable samples arise from two transmit PCM codes which when 
scaled by the same attenuation and quantized according to a thresholding procedure create identical PCM codes. 
An example taken from the attached attenuation table (pages A1-A2) which shows the attenuation for each uLaw 
25 PCM code when it is transformed by shifting from 1 to 26 codes numerically. 

Example: 



30 



35 



Transmitted 


Received shifted by 1 PCM code 


uLaw 


Linear 


uLaw 


Threshold 


Attenuation 


dB 


129 


7775 


129 


7647 


0.983537 


-0.14419 


128 


8031 




7903 


0.984062 


-0.13955 



Processing of Indistinguishable PCM Codes 

40 The processing of the feature set relies on knowledge of the coding law (uLaw or Alaw) used. In practice, it is nec- 
essary to consider the preciseness of the numerics used in the attenuation process. Suggestions for accommodating 
implementation deviations is discussed in a later section. The attached C language source code provided as part of this 
application (pages B1- B21) implements a demonstration of the identification procedure indistinguishably codes. The 
basic algorithm is the following: 

45 

1. For each indistinguishable PCM code received, the original pair of transmitted PCM codes can be determined 
since the PCM codes are sent in a known training sequence. A set of indistinguishable transmitted PCM code pairs 
are determined from this training sequence. 

2. For the first indistinguishable PCM code pair, determine the minimum attenuation required for the larger trans- 
50 mitted PCM code to be attenuated to a number of successive lower PCM codes. (A set of 30 successive lower PCM 

codes would be sufficient to cover an attenuation range of approximately 1 1 dB.) Then determine the maximum 
attenuation permitted for the smaller transmitted PCM code to be attenuated to the same successive lower PCM 
codes. (Minimum attenuation is numerically closer to one, i.e., a lesser attenuation, while a maximum attenuation 
is a greater attenuation.) If the maximum attenuation is a smaller attenuation than the minimum attenuation, discard 
55 the attenuation range as it represents a missing code rather than an indistinguishable code. The minimum and 
maximum attenuations for each successive lower PCM code forms an initial set of candidate attenuation ranges. 

3. For each successive indistinguishable transmit PCM code pair, another set of possible attenuation ranges is. 
determined using the procedure as described in Step 2. Each element of the candidate attenuation range set cre- 
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ated in Step 2 is searched for either full or partial overlap with any element of the possible attenuation range set. 
(In case of partial overlap, the attenuation range in the candidate set can be reduced for refined attenuation accu- 
racy.) In case of no overlap, the candidate attenuation range should be discarded. The discard occurs because two 
different indistinguishable codes must arise from the same linear transformation. 

4. Step 3 is repeated for successive pairs of indistinguishable transmit PCM code pairs until either a) the set of can- 
didate attenuation ranges is reduced to a single range or b) all the pairs of indistinguishable transmit PCM code 
pairs have been processed. 

In the first case, the actual attenuation is bounded by the single remaining attenuation range. We can say it is the 
detected attenuation can be the average of the high and low attenuations. Alternatively, if the history of each overlap- 
ping attenuation range is maintained (or regenerated) a probabilistic approach of determining a weighted median value 
for determining the detected attenuation can be implemented. 

For the second case where the attenuation appears to remains non-unique, the set of candidate attenuation ranges 
can be used to determine multiple sets of indistinguishable PCM code pairs can be determined for each candidate 
attenuation. The generated sets of indistinguishable PCM code pairs can then be matched with the original received 
set. The generated set using the correct attenuation should match nearly identically if not identically. The sets corre- 
sponding to incorrect potential attenuations would generate significant number of additional pairs of indistinguishable 
PCM code pairs or have missing indistinguishable PCM code pairs. Thus the correct detected attenuation can still be 
determined using this procedure. (Multiple candidate attenuations appear to arise from attenuations of 6 dB or greater 
combined with the decreasing step sizes used in with uLaw or Alaw companding.) 

Processing of Missing PCM Codes 

The procedure described above for using indistinguishable PCM codes to determine the attenuation can be 
adapted to operate with missing received PCM codes. The transmitted PCM codes corresponding each of the received 
PCM codes around the missing PCM code can be determined. Follow the same algorithm except that the minimum and 
maximum attenuations are exchanged in the algorithms. In other words, the larger transmitted PCM code produces the 
maximum attenuation and the smaller transmitted PCM code produces the minimum attenuation. With this adaptation, 
the algorithm described remains essentially the same. 

Example: 



Transmitted 


Recieved shifted by I PCM code 


uLaw 


Linear 


uLaw 


Threshold 


Attenuation 


dB 


144 


3999 


144 


3935 


0.983996 


-0.14013 


143 


4191 




4063 


0.969458 


-0.26942 



Other Variations 

For robustness with real signals corrupted by noise, it may be desirable to allow a certain number of pairs of indis- 
tinguishable PCM codes not to match a candidate range. An alternative to simply dropping a candidate attenuation 
range would be to use a probabilistic greatest likelihood model for selection of the detected attenuation. 

Furthermore, given that implementers of digital attenuation pads may understand and implement the mapping 
process slightly differently (i.e., use of G. 71 1 threshold values or average between PCM codes) and allow various inac- 
curacies in the numerical operations (numeric representation, rounding, truncation, etc.), each candidate attenuation 
range may need to be broadened to make overlap again more likely. The broadening can be by a fixed value or a rela- 
tive amount about the attenuation computed. The previous suggestion of discarding a certain number of pairs of indis- 
tinguishable PCM codes also addresses this problem as very few pairs are likely to be affected by these numerical 
inaccuracies. 

As for digital attenuation pads that pass through the robbed bit signaling (RBS) information, the algorithms pre- 
sented can also be applied except that v the indistinguishable transmit PCM code pairs will be separated by one code 
and a second indistinguishable transmit PCM code pair will adjacent. The candidate attenuation ranges become bigger, 
but still can be matched by the procedures presented in this invention. 

The procedures described in this invention suggest a robust manner to identify the digital channel attenuation 
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though the use of received PCM codes corresponding to a known transmit PCM code sequence. The use of such an 
algorithm would be necessary when the FCC permits the channel power to be increased in compensation for the sys- 
tematic digital attenuations. 

We further suggest that knowledge of the exact received PCM code set may be preferred by both the analog and 
5 digital PCM modems. Echo cancellation in the digital PCM modem may operate with incrementally better performance 
with the knowledge of actual PCM codes presented to the Codec. The receiver's desired spectral shaping can be more 
precisely honored by having the transmitter use the knowledge of the actual PCM codes presented to the Codec in the 
shaping function implementation, then reversely mapping the PCM code after shaping to the transmit PCM code to be 
sent. 

10 The exclusive use of pre-identrf ied attenuation table mappings would make the PCM modem technology incapable 
of operating with existing private telephone equipment. It is recommended that the attenuation pad mappings be dis- 
covered for each connection. 

Claims 

15 

1 . A method of determining digital channel attenuation.comprising: 

receiving a known training sequence of PCM codes, which PCM codes are subjected to said attenuation within 
said digital channel; 

20 quantizing said received known training sequence of PCM codes according to a predetermined thresholding 

procedure; 

identifying identical PCM codes created as a result of said thresholding procedure; and, 

determining said attenuation of the digital channel based upon said identification of identical PCM codes. 

25 2. A method of determining digital channel attenuation, comprising: 

receiving a known training sequence of PCM codes, which PCM codes are subjected to said attenuation within 
said digital channel; 

quantizing said received known training sequence of PCM codes according to a predetermined thresholding 
30 procedure; 

identifying PCM codes omitted as a result of said thresholding procedure; and, 

determining said attenuation of the digital channel based upon said identification of omitted PCM codes. 

3. A method of determining a digital channel PCM code transformation comprising: 

35 

receiving a known training sequence of PCM codes, which PCM codes are subjected to said PCM code trans- 
formation within said digital channel; 

quantizing said received known training sequence of PCM codes according to a predetermined thresholding 
procedure; and 

40 determining the transformation of transmitted codes to those received. 

4. A method of improved echo cancellation in a communication network having an analog and a digital modem.com- 
prising: 

45 saving codes transmitted from said digital modem to said analog modem for echo cancellation; 

transforming, by a mapping table, codes transmitted from said digital modem to codes received by said analog 
modem; and 

using said received codes as a reference signal for cancellation of echo. 

so 5. A method of improved spectral shaping using a transmit shaping transfer function in a communications network 
having an analog and a digital modem, comprising: 

transforming, by a mapping table, codes transmitted from said digital modem to codes received by said analog 
modem; 

ss using said received codes for transformation to their linear value equivalent representations; and, 

applying said linear value representations to said transmit shaping transfer function. 
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********* — - ****** */ 

pcmmap.c 

(C) 1997 VoCAL Technologies Ltd. 

ALL RIGHTS RESERVED. PROPRIETARY AND CONFIDENTIAL 

VoCAL Technologies Ltd. 
3032 Scott Blvd. 
Santa Clara. CA 95054 

Pioduct C 

Module. PCM 

This file contains (be PCM mapping functions. 

Revision Number: $Revision$ 

Revision Status: $State$ 

Last Modified $Date$ 

Identification. SldS 

Revision History: $Log$ 
Revision 1.0 1997/03/01 00:00 00 VD 
Initial release of software 



^include standard. h" 
^include "pcm.h" 
#include <stdio h> 
#tnctude <math.h> 



/* 



struct law_array_s { 
sint 15 code: 
sint 15 'inear; 
sint 15 other: 
sint 1 5 deselect: 

}- 

struct laws ( 

float tx_power_db: 
float rx power_dta: 
sint 1 5 count: 
sint ) 5 dmin: 

sint 15 indistinguishable: 
sint 15 x; 
sint3 1 bit rate. 

struct law array s value[128]: 
} ujaw. u law rx*" 

#defineLAW SELECTED 0 
#define LAW DESELECTED MIN 0x0001 
#define LAW DESELECT ED~MAX 0x0002 
tfdefine LAW DESELECTED DMIN 0x0004 
^define LAW'DESEt EC TED AVOID 0x0008 
#defineLAW DESELECTED DUPLICATE 0x0010 
^define LAW. DESELECTED_POWER 0x0020 

/ * * 

#define MAP FRAMF S1ZE 6 

uint48 frame_slot_modulus|MAP_FRAME_SIZE]: 
uint48 totat_symbols_per_frame, 

/ 

#define USE_4D_TRELLIS 



void 

u law in'it (struct law s *law) 
{" " 
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sint 15 idx, pcm; 

law->tx_power_db = 0.0; 
law->ix_power_db = 0.0: 
law->cou(it = 0; 
law->indistinguishable = 0: 
law->dmin = D. 
Iaw->bit_rate = OL; 

for (idx = 0: idx < 128; + + idx) { 
pcm - 255 - idx; 



law->value 
law- > value 
law->value 
law->value 



idx 
idx 
idx 

idx' 



code = pcm: 

linear = u law pcm decode (pcm), 
other = ulaw pcm threshold (pcm); 
deselect = LAW SELECTED; 



} 

r 



*/ 



#define SELECTLARGER 
void 

ujaw_attenuate (struct law s Maw. float attn_db) 



/* appears to work better with standard 3 and 6 dB attenuation V 



sint 15 idx, pern: 
float attn; 

attn = dbjojloat (-attn db); 

law->tx_power_db = 0.0: 
law->rx_power_db = 0 0; 
law->count = 0; 
law->indistinauishable = 0: 
law->dmin = "0; 
law->bit_rate - 0L; 

for (idx = 0; idx < 128. **-idx) { 
pcm = 255 - idx; 



other = (sint15) ((float) (u law_pcm_decode (pern)) * attn); 
.code = ujawpcmencocle (law->value[idx]. other); 
linear - u law pcm decode (law->value[idxj.code); 
deselect = LAW SELECTED; 



law->value[idx-1] deselect = LAW_DESELECTED_DUPLICATE: 



law->value[idx 
law->value idx 
law >value idx 
law->value[idx 
// if (idx »= 0) { 

// if (iaw->value[idxl code == law->value[idx-1].code) { 

//#ifdef SELECT LARGER 
// 

//#else 

law->value[idx|. deselect = LAW_DESELECTED_DUPLICATE; 



II 

//#endif 

II 
II 

) 



) 



} 



/* 



void 

u law range (struct law s Maw. sint 15 min, sint 15 max) 
sint 1 5 idx; 

for (idx = 0. idx < 128; + Hdx) { 

if flaw->value[idx). linear < min) { 

^ * law->value[idx). deselect i= LAWDESELECTEDMIN; 
if tlaw->value(idx] linear > max) { 

^ law->value[idx].dese!ect |= LAW_DESELECTED_MAX; 



*/ 



void 

u_law_dmin (struct law_s Maw. sint 15 dmin) 
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sint15 idx, pcm; 
sint15 current; 

current = law->value[125l. linear; 

law->value[127].deselect |= LAW DESELECTED_AVOID; 
law->value 126 deselect = LAWDESELECTED AVOID; 



for (idx = 124; idx >= 0; -idx) { 

if (law->value[idxl. linear > (current - dmin)) { 

law->value[idx]. deselect |= LAW_DESELECTED_DMIN; 



else { 

current = law->value[idx]. linear; 

void 

ujawdiscardjndinstinguishable (struct law_s *law) 
sint15 idx; 

for (idx = 1; idx < 125; ++idx) { 

if (law->value[idx].coae == law->value[idx-1].code) { 
//#ifdef SELECT LARGER 

// ~ law->value[idx-1J.deselect |= LAW DESELECTED DUPLICATE; 

//#else ~" ~~ 

// law->value(idx].deselect |= LAW DESELECTED DUPLICATE; 

//tfendif 

if (law->va!ue[idx-1].deselect & (LAW DESELECTED DMIN | 
LAW DESELECTED AVOID | LAW DESELECTED POWER)) { 

law->value[idx-1J.dese!eci |= LAWDESELECTE DO U PLICATE; 

else { 

law->value[idx].deselect |= LAW_DESELECTED_DUPLICATE; 

, • ' ' 

I* ****-*********************************#*********************************** * i 
void 

ujaw_attenuate_dmin (struct law_s *txjaw, struct !aw_s *rx_law) 
sint15 idx; 

rx law->value[127].deselect|=LAW DESELECTED AVOID; 
rx~law->value[126].deselect |= LAW^DESELECTED~AVOID; 



for (idx= 124;idx>=0; -idxW 
-- r " *- v->valuefidx].d 

rx_law->value[idxl.deselect |=H_AW_DESELECTEDJ 



if (tx law->valuef idx] deselect & LAW DESELECTED DMIN) I 

DMIN; 



******* */ 



, • 1 

7* A*****************************************************************: 

void 

u Jaw exclude (struct laws *law) 

sint15 idx, pcm; 

for (idx = 127; idx >= 0; --idx) { 
pcm = 255 - idx; 

if ((pcm == 165) || (pcm == 169)) { 

law->valuelidx].deseiect |= LAWDESELECTEDAVOID; 

• 1 ' 

I* ************************************************************************* * i 
float 

u_law_power (struct law_s *law, float max_db) 
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sint 1 5 idx. count, dmin, prev; 

float factor, power, rnax_float; 

float tx_power tx_sum, rx_power, rx_sum; 

factor = pow (10.0, ((3.17 + 3.01244) / 20.0)) / 8159.0; 

count = 0; 
tx_sum = 0.0; 
rx_sum = 0.0; 
dmin = 8159: 
prev = -8159; 

max float = db_to_power (max_db): 
for (idx = 0; idx < f28; ++idx) { 

if (law->value[idx]. deselect == 0) { 

rxjDower = law->value[tdx]. linear * factor; 

rx_power = rx_power * rx_power; 

rx_power = rxpower + rx_sum; 

tx_power = ujaw.value[idx]. linear * factor; 
tx_power - tx_power * tx_power; 
tx_power = tx_power + tx_sum; 

power = tx_power / (float) (count + 1); 
if (power <= max_float) { 

rx_sum = rx_power; 

tx_sum = tx_power; 

count++; 

if ((law->value[idx]. linear - prev) < dmin) { 

dmin = law->value[idx]. linear - prev; 

prev = law->value[idx]. linear; 

else { 

^ law->value[idx].deselect = LAW_DESELECTED__POWER; 

//printf ("%f %f Vod %d\n'\ power, tx_sum. count, law->value[idx].deselect); 

ix_power = tx_sum / (float) count; 
tx_power - power_to_db (tx_power); 

rx_power = rx_sum / (float) count; 
rx_power = power_to_db (rx_power); 

law->tx_power_db = tx_power; 
law->rx_power_db = rx_power; 
law->count = count; 
law->dmin = dmin; 

return power; 

y* ******************************************************************** *f 

void 

u_law_count_indinstinguishable (struct law_s *law) 
sint15 idx, count; 
count = 0; 

for (idx = 1; idx < 128; ++idx){ , 
if <law->value[tdx]. deselect == LAW_DESELECTED_DUPLICATE) { 
count++; 

iaw->indistinguishable = count; 

} 

t* ************************************************************************* *j 
sint31 

ujaw_bit_rate (struct law_s *law) 
^ float bit rate float; 



bit_rate float = (log 10 ((float) law->count * 2.0) I Iog10(2.0)) * 8000.0; 
law->biCrate = (sint31) bit_rate_float; 
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return law->bit rate; 



************* 



//tfdefine MAX SHIFT 20 /* covers to -7 08 dB 7 

#define MAX SHIFT 30 /* covers to -1 1 2 dB V 

//#define MAX_SHIFT 40 /* covers to A A 65 dB V 

float 

uJaw_detect_attenuation (struct law_s law) 

sin(15 idx. idx1. idx2. idx3. trim; 
sint15 count: 

float attn_low[MAX SHIFT], attn high[MAX SHIFT! 

oat a lowfMAX_SHIFTJ. a highTMAX SHIFT]; 
float afln. cfiff; ~ ~ J 

count = 0; 
idx3 = 0; 



for (idx = 128; idx >= MAX SHIFT' -idx) { 
if (law->value[idx|. deselect & LAW 
idx3++; 

if (count -= 0) { 



DESELECTED DUPLICATE) { 



#ifdef SELECT_LARGER 

u_law.value[idx+1). linear; 

u_law.value[idx]. linear; 
#e!se 

u_law.value(idx|. linear; 

u law. valuefidx- 11. linear: 
#endif 
// 

attn_high[count], attn_low[count^) 



for (count = 0; count < MAX SHIFT ++count) { 

attn_high[count] = (float) uJaw.value[idx-count], other / (float) 
attnJow[ count] = (float) uJaw.value[idx-counM [other / (float) 



attn_high[count] = (float) uJaw.value[idx-count-1].other / (float) 
attn_low[count] = (float) uJaw.value[idx-count-2J.other / (float) 



} 

else { 

#ifdef SELECT_LARGER 

ujaw valuefidx + 1]. linear; 

ujaw.valuefidxj. linear; 
#else 

u_law.value[idxj. linear; 

u law.value[idx-1l. linear 

#endif 

// 

ajiighpdxl], ajowfidxl]); 



printf ("original %d %f %f\n", u Jaw. valuefidx]. code, 
count = MAX_SHIFT; 



a_high[idx2])) { 
attn_high[idx1], attn_low[idx1], a_high[idx2]) 



aJow[idx2])) { 

attn_high[idx1J, attnJow[idx1], a_low[idx2]); 



for (idxl = 0; idxl < MAX_SHIFT; ++idx1) { 

a_high[idx1] = (float) u Jaw. value[idx-idx1 ]. other / (float) 
a_low[idx1] = (float) ujaw.valuefidx-idxl -1J. other / (float) 

a_high[idx1] = (float) ujaw.value[idx-idx1-1|. other / (float) 
a_low|idx1] = (float) u_law.value[idx-rdx1-2].other / (float) 

printf ("generating %d %f %An", u_iaw.value{idx].code. 

} 

for (idxl = 0; idxl <= count; ++rdx1) { 
trim = 1; 

for (idx2 = 0; idx2 < MAX SHIFT; +-Hdx2) { 

if <(aUn_high[idxfJ>= a_high[idx2]) && (attn_low[idx1] 



printf ("reduce high %f %f %f\n M , 



attn high[idx1] 
1rim~=0; 



: a_high[idx2]; 



if ((attn_high[idx1) >= aJow[idx2]) && (attnJow[idx1J <= 
printf ("increase low %f %f %f\n M , 



} 



attnJow[idx1] = a_lowlidx2]; 
trim = 0; 
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if ((atUi_!Aig:i[idx1] <= o high[iGx21) £ ■& (rttn IcMidxt] 
a Jow[idx2])) { * 

printf ("hold r ,m %f\n\ attnjngn[idx I J. 



attnjowjidxl]): 



, > ' 

if (count <- 0) { 

printf ( count break Vn"): 
break; 

\\ (idx3 >= 10) break: 



trim = 0; 

(attn_high[idx 1] < attn_low(idx1]) { 
trim = 1 ; 

if (trim) { 

printf ("trim %f %f\n". attnjitghfidx 1], attn low[idx1]); 

for fidx2 = idxl ; idx2 < (count - V); + * idx2) {" 
attn_high[idx2] = attn highlidx2 + i]: 
attnJow[idx2] = attnjbw[idx2 H]: 

-count: 
-idxl: 



\ (count > 1) { 



attn = (attn higlifO] + attn lowfOJ) / 2.0: 
diff - attnjfigh[d) - aitnJow[0J: 
for (idxl = 1 : idxl < count, + + idxt) { 

if «attn_high(idx1] - attn low[idx1J) < diff) { 

attn = (attn highlidxi) * attnjowjidxl]) / 2.0. 
diff = atlnj?igh[idx1] - attn_low(icfx1 ]: 

//v printf f Muplitpte potential attenuations detected, resolve by generating and matchingW 

//v "duplicated codes to those observed. Currently choosing largest 

range. \n"); 

else { 

attn = (attn_high[0] + attnJowfO]) / 2.0: 

attn = loglO(attn) * 20.0: 
// printf ("count %d attempts %d range %f v 0 f\rf count. idx3. attnjiigh[0). attnjow[0]); 

// printf ("attenuation detected 3 /of\n". attn): 

return attn. 

} 

r — * — * / 

sint 15 

u Jaw_map. symbol set (struct taw_s 'Ixjaw. struct law s 'rxjaw) 

sint 15 idx. count, limit: 
uintIG mask. 

total symbols _pei_fiame.lsw = 1. 
totarsymbols pei_ frame. mid = 0; 
totals symbols" per frame. msw - 0. 

// printf("%0<1x %0^1x %0 / lx\n ,, [ total_symbols per frame. msw, 

// totat_symbols_perJiame.mid. toTat_symbols_per Jrame.lsw); 

for (idx = 0; idx < MAP__FRAME SIZE; +-+idx) { 



frame_slot_modulus 
frame_slot_modulus 
frame slot modulus 



idx 
idx 
idx 



Isw = total symbols_per Jrame.lsw; 
mid - totarsymbols_perjrame.mid; 
msw = totaT_symbols_pe7jrame.msw; 



#ifdef USE_4D TRELLIS 

"limit = tx law->count; 
if (idx & T) { 

limit = limit » 1 ; 
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for (count = 1; count < limit; + + count) { 

uint48_add (&total_symbols per_frame. &frame_slol_modulus[idx]. 
&tofal_symbols_per_frame); 

#else ^ 

for (count = 1; count < tx law->count; ++count) { 

uint48_ add (&totaf_symbols per_frame, &frame_s(ot_modulus{idx]. 

&tofal_symbols_per_ftame); 

#endif * 

// printf("%04x %04x %04x\n". tota1_symbols per_Jrame.msw. 

II ^ total_symbols_per_frame.'mid. toFal_symbols_per_frame Isw): 

count = 47: 

mask = 0x8000: 

for (idx = 0: idx < 16; ** + idx) { 

if (mask & total_symbo!s per frame msw) { 
return count; 

mask = mask >> 1: 
-count; 

} 

mask = 0x8000; 

for (idx = 0; idx < 16: +-Hdx) { 

if (mask & lotat_symbols_per frame. mid) { 
return count: 

mask = mask >> 1 ; 
--count; 

} 



mask = 0x8000; 

for (idx = 0; idx < 16; + + idx) I 

if (mask & tatalsyinbols per frame. Isvv) { 



return count; 

mask ~ mask >> 1 : 
--count; 

} 

return count: 



void 

Drint deselect (sint 15 deselect) 

if (deselect & LAW DESELECTED_MIN) { 
^ printf ("(Below M;n) *'); 

if (deselect & LAW_DESELECTED_MAX) { 
printf ("(Above Max) "): 

if (deselect & LAW_0ESELECTED_AVOID) { 
^ printf ("(Avoid Dmin) "); 

if (deselect & LAW DESELFC1 EDDMIN) { 
printf ("(Refween Dmin) "): 

} 

if (deselect & LAW_DE SELECTED DUPLICATE) { 
printf ("(Duplicated Code) "Y. 

if (deselect & LAW_DESELEC T ED POWER) { 
^ printf ('(Power Limit) "); 

^ printf("\n"); 
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//#define DO_FEW 

/ ***** - */ 

//tfdefine DO_PRINT 0 
//^define DO_PRINT*~1 
//^define DO_PRINT~2 
//#defrne DO PRINT"3 
//#define DO"PR!NT~4 
//#define DO PRINT~5 



#ifdefDO FEW 
#define MAXATT N 1 
^define MAXDMIN 4 

^define START DMIN 48 
#define STEP DMIN 4 

tfdefine START At IN M O 
#define STEP ATTN < 0 



r 



#else/* DO FEW V 
^define MAX ATTN 14 
#define MAX DMIN 36 

tfdcfinc START DMIN 4 
fldefine SI EPDMIN 4 

tfdefine START ATTN 0 0 
#defineSTEP ATTN 1.0 
#endif r DO FEW 7 



/' 



void 

main (void) 
{ 

smt15 idx, count; 
float attn. 

sint 1 5 bits per_frame; 

sint31 bit rate. overal1_bit rate: 

float input attn, 

sint 1fi input txdmin. input rx dniin. 
sint 15 i. j. results idx. 

struct results_s { 

float attn; 

sint 15 dmin. 

sint31 bit_rate_simple; 

sint3l bit rale ^exclude; 

sint31 bit rate ~distinguishable_tx_drnin. 

sint31 bit rate distinguishable rx dmin. 
} -esults [MAX ATTN T MAXiDM(N); 

resultsjdx - 0; 

;nput_attn = START ATTN; 
// input Jxdmin = 1 2CTT 
// mput_rx_dmin = 120; 

for (i = 0; i < MAX ATTN, input attn += STEP ATTN. i + + i { 
: nput_tx_cTmin = S I ART DMIN; 
input rx dmin = START" DM IN; 

for (j = 0: j < MAX^DMIN, input Jx__dmin += STEP_DMIN. input_rx_dmin += STEP_DMIN. j++) { 

printf ("\n"). 

u taw_init (&ujaw): 



// ujaw_jange [&u law. 30. 4000) 

u lawdmtn (&u_Taw. input tx dmin); 

// uJaw_power {iujaw, -12.D /^9.479817); 
u_law_power (Sujaw, -12.0); 
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ujaw_bit_rate (&u_lavv); 

prinlf ("%6.2f %3d %3d ". input__attn. inputjx__dmin. input_rx_dmin); 

printf ("tx power %3.2f rx power %3.2f count %d dmin %d dup %d bit rate %ld\n". 
prinlf ("tx pwr %3.2f rx pwr %3.2f cut %2d dmin %3d dup %2d rate %ld\n". 

ujaw.tx powered b. u_law.rx_power db, ujaw. count, ujaw dmin. 

u law. indistinguishable. ujaw. bit rale); 



*ifdef DO_PRINT 0 

for (idx = 0. idx < 128; ++idxH 

" 5d %5d %5d %5d '. u_ law. valuefidx). code, u law.value[idx].linear. 
ujaw. value[idx]. other, ujaw pern encode(u Taw valuejidxj.lineai ). 



tfendif 



} 



u law pem encodefujaw.valuefidx], other - f)); 
print deselect (Lf law. vafuejidx]. deselect): 



u law_attenuate (&ujaw rx. input attn); 
A' ujaw_range (&u Jaw_rx.~22. 40007: 

// ujaw_dmin (&u law_rx. input_rx dmin); 

Lj_law_atteruiate"dmin (&u law. Slijaw rx) /* transfer from transmitter 7 
ft ujaw_exclude (&ujaw_rx^ 
// u_law_discard tndinstinguisliable (&u law_rx). 

uJaw_power f&u law rx. -12.0); 

ujaw_discardjn3insfihguishable (&u_law_rx) 

u lawcounl indistinguishable (&ujaw_rx). 

ujaw_bit rate (&u law_rx); 

printf ("%G 2f %3d %3d ". input_attn. input tx dmin. input_rx dmin); 

// printf ("tx power %3 2f rx power %3.2f count %ti dmin %d dup %d bit rate %ld\n". 

printf ("tx pwr %3.2f rx pwr %3.2f cnt %2d dmin %3d dup B / C 2d rate %ld\n". 

ujaw_rx tx_power_db. ujaw rx rx power_db. ujaw rx count, 
ujaw rx.dmin, u law_rx. indistinguishable, ulaw rx.bTt rate); 

#ifdef DO_PRINT 1 

for {idx = 0; idx < 128; ++idxW 

printf ("%3d %5d %5d %5d ", u law.value[idx].code, ujaw rx.value[idx]. other. 

u law_rx.value[idx]. linear. uJaw_rx.value[idx7code); 
print_deselect(u_Iaw_rx value[idx] deselect); 



flendif 



attn = u law_detect attenuation (&u iaw rx). 
printf ("attenuation detected %f\n". altn): 



resultsfresultsjdx] attn "input attn; 
results[resultsjdx].dmin = inpu7jxjjmin: 
if (u_law_rx dmin == 0) { 

results[results idx] bit_/ate__simple = 0; 

else { 

^ resiilts[results_idx].bit_rate_simple = u_law_rx.bit_rate; 



u law allemiate (&ulaw rx. input attn); 
// u, lawjange (&u law_rx.~22. 40007; 

// ujaw_dmin (&ujaw_rx, input_rx dmin); 

LJjaw attenuate dmin (&u law. Sujawrx); /* transfer from transmitter 7 

u Iaw_exclude (5u law_rxj; 
// u_law_discard indistinguishable (&u_law_rx): 

u law power (&u law rx. -12.0); 

u Jaw^discard Ji iclinstu iguisha ble (&u_law_r x) . 

uJaw_count indistinguishable (&ujaw_rx); 

u_law_bit_rale (&u_law_rx); 

printf ("%6.2f %3d %3d '\ inputattn. inputjx_dmin. input_rx_dmin); 

// printf ("tx power %3.2f rx power %3.2f count %d dmin %d dup %d bit rate %ld\n'\ 
printf ("tx pwr %3.2f rx pwr %3.2ff cnt %2d dmin %3d dup %2d rate %ld\n". 



u law_rx.tx power_db. u law rx.rx power_db, ujaw rx. count, 
ujaw_rx.dmtn. ujawrx Indistinguishable. ujaw_rx.r5it_rate); 



#ifdefDO PRINT_2 
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for <idx = 0: idx < 128; ++idx) { 
%5d 



} 



printf ("%3d %5d %5d %5d u law. valuefidx]. code, ujaw rx.value[idxj. other. 

u law_rx.value[idx). linear, uJaw_rx.value[idxJcode); 
print_dese!ect(u_faw_rx.\/alue[idx). deselect); 



#endif 



at(n = u law_detect attenuation (&u faw_rx); 
printf ("attenuation detected %f\n". altn); 

if fujaw_rx.dmin == 0) { 

resulls[resultsjdx|.btt_rate_exc1ude = 0: 

else { 

results[resiMsjdx).bit_rate_exclude = u law_rx.bit_rate; 



ti law_attenuate (&u_law rx. input attn); 
// u law_rnngc (&u law_rx.~22. 4000); 

// u_law_dmin (&ujaw_rx. input_rx dmin); 

ujaw_attenuate_dnun (&ujaw. £1/ law_rx). 

u_law_discard indinstinguishable (Sulawrx): 

ujaw power fS.tiJaw_.rx. -12.0): 
// u_law_power (&ujaw_rx, (-12 0 - attn)); 
// uJaw_counl indinstinguishable (&ujaw_rx); 

ujaw_bit_rafc (&ujaw_rx); 

printf C'%6 2f %3d c/ n3d M . /nput_attn. input Jx_dmin. input_rx_dniin); 

// prinif ("tx power 0/ c 3 2f rx power %3 2f count %d dmin %d dup %d bit rate %Id\n" 
printf ("tx pwr %3.2f rx pwr %3 2f cnt %2d dmin %3d dup %2d rate %!d\rV. 

n law_rx.tx_power_db. ujaw rx rx_power_db, ujaw_jx count 
u_law_rx.dmin. uJawj-xmdisTinguisnable. ujaw rx.bit rate); 

#ifrief DO_PRfNT 3 

for (idx = 0: idx < 1 28; ++idxl { 

printf ("%3d %5d %5d %5d u law value[idx]. code, ujaw rx.value[idx].other, 

u !aw_rx.value[idx). linear, ujaw rx.value[(dxf code); 
print_desetect(u_Taw_rx.value[iax]. deselect); 



#endif 



attn = u (aw detect attenuation (&u law_rx): 
printf ("altenuation detected %f\n". altn); 

results|resiiltsjdx] oit_rate_distinguishab!eJx_dmin = ujaw_rx.bit_rale. 
/ * ****** * 7 

ujaw_attenuate (&u_law rx. input attn); 
// ujaw_rancje <&u lawjx~22. 40007; 

ulawdmin <&ujbw rx. inputjx dmin); 
// u_law_attenuale_dmih (&ujaw. STu law_rx); 

u Jawdiscard indinstinguishable (8TuJaw_rx): 

u_law_powei ['&ujaw_rx. -12.0): 
// u law power (&ujaw_rx. (-12.0 attn)); 
// u_law"count indinstinguishable (&ujaw_rx) 

u Jaw "hit ..jafe (8u lawrx); 

printf <"%6 2f c «3d %3d inputattn, input Jxdmin, inputrxdmin); 

// printf ("tx power %3.2f rx power %3.2f count %d dmin %d dup %d bit rate %ld\n". 

printf ftx pwr %3 2f rx pwr %3.2f cnt %2d dmin %3d dup %2d rate %ld\n'\ 

l U aw _ rx tx_power_db. u law rx.rx Dower db, ujaw_rx. count. 
ujaw_rx.dmm. ujaw rx.TndisTinguisnable. ujaw_rx.bit_rate); 

#ifdef DO_PRINr 4 

for (idx - 0: idx < 128, + + idxi { 

printf r%3d %5d %5d %5d M . u law. vaiuefidx]. code, ujaw rx.value(idx]. other, 
u !aw_rx.value(idx]. linear, u Jaw__rx.value[idxJxode); 
^ print_deselect(u_faw_rx.value[idx], deselect); 

attn = u faw_detect attenuation (&u fawjx); 
printf ("altenuation detected %f\n", altn); 

#endif 

B10 
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resul!s(resutts_idx++] bit_rate_distinguishab!e_rx_dmin = ujaw_rx,bit_rate: 

/ * — — " / 

bits_per_fmme - uJaw_map_symbo!_set {&u_law_rx, &ujaw_rx): 

bit_rate = (sint31 ) (((float) (bits_per_frame) / 6.0 + 1.0) * 8000.0): 

#ifdef USE_4D TRELLIS 

overall^bit^rate = (sin*31) (((float) (bits_per_frame) /6.0+10 + 0.5) * 8000.0) 

#else 

overall_bit_rate = (sin:31) (((float) <bits_per_frame) / 6.0 + 1 .0) * 8000 0): 

#endif 

ffifdef DO_PRINT_5 

pnnlf ("bits per frame %6 bit rate %ld overall bit rate %ld (with trellis)\n'\ bits_per_fr3mc. 
bit rate, overall bit rate); 
#endif " ~ " 



r **" 
#if 0 



#endif 



for (coirtt = 0; count < 1 000; count = count + bits_per_frame) { 
generate_bits (bits per_frame): 
ujaw_send_frame (&u law. &u law__rx): 
u law_receive_fraine (STuJaw. £ujaw_rx): 
cfTeck_btts (bits_per_frame). 



r 



tfifndefDO TEW 



#endif * 



fonidx = 0; idx < MAX ATTN * MAX DMIN: idx++) { 

printf("%6.2f %3d %5ld %5lcP/o5ld %5ld\n' . resultsfidxj.attn, resultsfidxj.dmin. 

"-' ;j "Vbit_rate_simple. results[idx).bit rate_excluae. 

bit_rate_disttnguishable_tx_dnTin, 
. bit_rate_distingu isha blerxdmin); 



results 
results 
results 



idx 
idx 
idx 
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r 



******************* 



********* * J 



pcmcnv.c 

(C) 1997 VoCAL Technologies Ltd. 

ALL RIGHTS RESERVED. PROPRIETARY AND CONFIDENTIAL. 

VoCAL Technologies Ltd. 
3032 Scott Blvd. 
Santa Clara, CA 95054 

Product: C 

Module: PCM 

This file contains the PCM conversion functions. 



Revision Number: 
Revision Status 
Last Modified: 
Identification: 



$RevisionS 
$State$ 
$Date$ 
Sld$ 



Revision History: $Log$ 
Revision 1.0 1997/03/01 00:00:00 VD 
Initial release of software 



*/ 



********* 



*********** 



******* ****** ****** 



* 7 



^include "standard h" 
#include "pcm.h" 



/* 



#include <stdio.h> 
sint15 

ujawpcmencode (sfract15 input) 
sfract15 value, seg: 
value = input; 

if (input < 0) value = -input; 

value = value + 33; 

if (value > 8159) value = 8159; 

seg = iexp (value); 
;/printf( M %5d %5d %5d". input, value, seg); 
value = norm (value, seg); 

value = value A 0x4000; 
value = value » 10; 
//printf(" %5d". value); 

if (input < 0) { 

value = value I 0x0080; 

} 

seg = seg + 9; 
seg = seg « 4; 

value = value I seg; 
//printff %04x'\ value): 

value = value A OxOOff; 
//printf(" %5d\n", value); 
^ return (sint 15) value; 



#ifdef DO_U_LAW 

x 

x 

X 
X 
X 
X 
X 



AR = ABS AR; 

AY0 = 132; 

AR = AR + AY0; 

SR1 = 32636; 

IF AVAR = PASS SR1 

SE = EXPAR (HI); 



/* Take absolute value */ 

/* Limit values */ 

/* Compute exponent */ 

I* Normalize value */ 

/* Clear bit V 
/* Position bits */ 

/* Insert sign bit */ 

/* Compute segment */ 
/* Insert segment bits */ 
r Invert bits 7 



/* Input sample passed in AR (1 . 1 5) 7 
/* Take absolute value 7 
/* 33 « 2 7 

/* Limit values 7 

/* Compute exponent 7 
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x 

X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

#endif 



r Normalize value 7 



r Clear bit 7 
I" Position bits 

/* Create sign bit 7 



SR = NORM AR(LO); 
AYO = 0x4000; 
AR = SRO XOR AYO; 
SR = LSHITT AR BY -10 (LO): 
AXO = SF. AR = PASS AYO. 
ir POS AR = PASS 0. 
SR = SR OR LSHIFT AR BY -7 (LO) 
AYO = 7; 

AR = AXO + AYO. 

SR = SR OR LSHIFT AR BY 4 (LO): 
AYO = OxFF; 
AR = SRO XOR AYO: /* Invert bits 7 

RTS. r Pass on rate converted value 



I* Insert sign bit V 
/* 9 if using 8031 range 7 
/* Compute segment 7 
r Insert segment bits V 



r 



sfract 15 

u_law_pcm_decode (uint8 code) 

sfract 15 value, seg: 

value = code & OxOOff: 
value = value A OxOOff; 

seg - value & 0x0070: 
value = value A seg: 
seg = seg >> 4; 

if ((value - 0x0080) >= 0) { 

value = value 0xff80 
value = value + value + 
value = value « seg; 
value = 33 - value: 



else { 



} 



T Mask bits 7 
T Invert bits 7 

/* Isolate segment bits 7 
/* Remove segment bits 7 
r Shift to tow bits 7 

/* Determine sign 7 

r Process negative value 7 
33: r Add segment offset 7 

/* Position bits 7 

r Remove seament offset 7 



value = value + value + 
value = value << seg: 
value = value - 33; 



33: 



/* Process positive va'ue 7 
r Add segment offset 7 
/* Position bits 7 

r Remove seament offset 7 



return value. 



r 



sfract15 

u Jaw_pcm_threshold (uint8 code) 



} 



sfract 15 value, seg. 

value = code & OxOOff; 
value = value A OxOOff. 

seg = value & 0x0070: 
vafue = value & OxOOOf; 
seg = seg >> 4; 
seg+ * ; 

value = value * 1 7; 
value = value « seg. 
value = value - 33; 

if ((code & 0x80) == 0) { 
value = -value: 

} 

return value. 



#ifdef DO__A_.LAW 
aJaw_pcm_encode: 



/* Mask bits 7 
r Invert bits 7 

r Isolate segment bits */ 
/* Remove segment bits 7 
/• Shift to low bits 7 



/* Add segment offset V 
/' Position bits 7 

I" Remove segment offset 7 

/* Determine sign 7 

/' Process negative value 7 



AR = ABS AR; 
AYO = 511; /M27 7 
AF = AR-AYO; 
IF GT JUMP ajaw_pcm_enc_1 



/* Input sample passed in AR (1.15) 
r Take absolute value 7 
/* Check for zero segment 7 
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a_law_pcm_enc_ 1 : 



SR = LSHIFT AR BY -4 (LO) 
AR = 0x4000: 
IF NEG AR = PASSO: 
SR = SR OR LSHIFT AR BY -7 fLO); 
AYO = 0x55; 
AR = SR0 XOR AYO: 
RTS: 



/* Downshift bits 7 



/* Create sign bit 7 
" /* Insert sign bit */ 



/* Invert bits 7 

/* Pass on rate converted value 7 



SE = EXP AR (HI). r Compute exponent 7 

AXO = SE. SR = NORM AR (LOl: /* Normalize value V 

AYO = 0x4000; 

AR = SRO XOR AYO: /* Clear bit 7 

SR - LSHIFT AR BY -10 (LO): /* Position bits 7 
AR = PASS AYO: 

IF NEG AR = PASS 0: /* Create sign bit 7 

SR = SR OR LSHIFT AR BY -7 (LO); /* Insert siqn bit 7 

AYO = 7; 

AR = AXO h AYO; /* Compute segment 7 

ir I I AR = PASS 0: 

SR = SR OR LSHIFT AR BY 4 (LO); /* Insert segment bits */ 
AYO = 0x55: 

AR = SRO XOR AYO: /* Invert bits 7 

RTS: /* Pass on rate converted value 7 



ident(IDENT RX B1_ PCM DECODE): 



a_law_pcm_ decode: 



/* Input sample passed in AR 7 

/* MR^k bits for inversion 7 
r Invert bits 7 

r Set sign bit 7 
r Set LSB of interval 7 
12 (LO): /* Isolate segment and interval */ 



#endif r DO A LAW 7 



AYO = Dxftt 
AR = AR XOR AYO; 
SR1 = 0x8; 
SRO = 0x800: 

SR = SROR LSHIFT AR BY 
AYO = 32; 

AXO = AR. AF = PASS AYO: 
AYO = 9; 
AR = SRI - AYO; 
IF LT AF = PASS 0; 
IF EQ AR = PASS 0; 
SR = LSHIFT SRO BY -11 (LO): 

SE = AR. AR = SRO OR AF; /* Add bit if necessary 7 
SR = LSHIFT AR (LO); r Position output 7 
SR - LSHIFT SRO BY 3 (LO): 
AYO = 0xFF80: 

AR = SRO. AF = AXO + AYO; r Check if sign b ; t set 7 
IF LT AR = - SRO; /* If it is. negate result 7 

RTS: /* Pass on rate converted value (1. 



/* Segment bias 7 
I* Determine shift 7 
r No extra MSB bit 7 
/* No less than zero bits 7 
Isolate interval 7 



15) 7 
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r * * - — — v 

r 

malh.c 

* (C) 1997 VoCAL Technologies Ltd. 

; ALL RIGHTS RESERVED. PROPRIETARY AND CONFIDENTIAL 

VoCAL Technologies Ltd. 
3032 Scott Blvd 

* Santa Clara. CA 95054 

\ Product. C 

I Module. MATH 

* This file contains fract ; onal math support functions 

Revision Number $Revision$ 

Revision Status* $State$ 

Last Modified $Date$ 

* Identification: $ld$ 

Revision \ listory. $Log$ 

Revision 1.0 1997/03/01 00.00:00 VD 

Initial re ! ease of software 

7 

/ * * * V 

^include "standard h" 
^include "pcm.h" 
^include <rr:a!h.h> 

/* ~ * * *** 7 

sint15 

iexp (sft act 1 5 value) 

sint 1 5 exp. 

if (value < 0) { 

value = -value: 

} 

exp- 15; 

if (value == 0) { 
return exp: 

} 

exp = 0: 

while ((value & 0x4000) == 0) { 

--exp: 

value = value « 1 ; 

) 

return exp 



sint 15 

norm (sfract15 value, sint 15 exp) 
return (value « exp): 



r 



float 

dbjojloat (float db) 

* float db_float: 

db float = db / 20.0; 

dbjloat = pow (10.0. dbjloat); 

return db float; 
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float 

dbjo power (float db) 

' float dbjloat; 

dbjloat = db / 10.0; 

dbjloat = pow (10 0. db float); 



return db float: 



float 

float_to_db (float power) 
* float float db. 



} 

r 



float jib = Iog10 (power) • 20 0: 
return float db: 



float 

oowertodb (float power) 
float power_db; 



} 

r 



power db = log 10 (power) ' 10 0; 
return power db; 



void 

uint48_add (uint48 *src1. uinl48 *src2. uint48 *dst) 

uint16 s1f slm. slh. s2l s2m. s2h. dstl. dstrn. dsth. 

s 11 = src1->lsw; 
slm = srd ->mid. 
slh - srd >msw; 

s2l = src2->lsw; 
s2m = src2 >mid; 
s2h = src2->msw: 

asm { 

rnov ax. s1l 
add ax. s2l; 
rnov dstl. ax; 

mov ax. s 1 m; 
adc ax. s2tn. 
mov dstm. ax; 

mov ax. slh: 
adc ax. s2h; 
mov dsth. ax. 

} 

dst->lsw = dstl: 
dst->mid = dstm. 
dst->msw = dsth; 

// printf("%04x %04x %04x = %04x %04x %04x + %04x %04x 0 /c04x\n" 

// dsth. dstm. dstl. slh. slm. s1l. s2h. s2m. s2l); 

) 

r — — v 

void 
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umt48_sub (uint48 *srd. uin!48 *src2. uint48 *dst) 

uint16 s1l. s1m. s1h. s2l. s2m. s2h. dstl. dstm. dsth: 

s1 1 = src1->lsw; 
s1 m = src1->mid. 
slh = srcl->msw. 

s2l = src2->lsw; 
s2m = src2->mid. 
s2h = src2->msw; 

asm { 

mov ax sll: 
sub ax. s2l; 
mov dstl. ax; 

mov ax. si m: 
sbb ax. s2m; 
mov dstm. ax; 

mov ax. s1h: 
sbb ax. s2h; 
mov dslh. ax: 

} 

dst->lsw = dstl: 
dst->mid = dstm: 
dst->msw = dsth: 

> 



sfract15 

frnpy (sfractl 5 a. sfractl 5 b) 

long a long, b long; 
short a_short: 

a long ~ a; 
b_long - b: 

a long = a long * bjong: 
ashort - (short) (a_long » 15); 

return a short: 
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— — * — — — — ^ 

pcm.h 

(C) 1997 VoCAL Technologies Ltd. 

ALL RIGHTS RESERVED. PROPRIETARY AND CONFIDENTIAL. 

VoCAL Technologies Ltd 

3032 Scott Blvd 

Santa Clara. CA 95054 

Product: C 

Module: PCM 

This file defines the PCM utility functions. 



Revision Number 
Revision Status 
Last Modified. 
Identification: 

Revision History: 



$Revision$ 
$State$ 
$Date$ 
$ld$ 

$Log$ 



#ifndef 
^define 



PCM PCM H 
PCM PCM M 



/* ***** *«. 

typedef struct { 

uint16 Isw, 
uint16 mid; 
uint16 msw: 

} uint48: 

sint15 iexp (sfract15 value). 

stnt15 norm (sfract15 value, sint 1 5 exp): 

float dbjojloat (float db): 

float db Jo power (float ab): 

float float Jo db (float power): 

float power Jb_db (float power): 



void uint48 



add (i 

rsub (t 



iin!48 *srd 



uint48 *src2. 
uint48 *src2. 



uint48 
uint48 



dst); 
dst): 



sintlS u lav/_pcm_encode (sfract15 input): 
sfract15~u law pcm decode (uintS code); 
sfract15 u law pcm threshold (uint8 code): 



r 



#endif /' PCM PCM H V 
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/* •* — — ****** — **** v 

* standard. h 

* (C) 1994. 1995. 1996 VoCAL Technologies Ltd. 

; ALL RIG! ITS RESFRVFD. PROPRIETARY AND CONFIDENTIAL. 

VoCAL Technologies Ltd. 
3032 Scott Dlvd 
; Santa Clara. CA 95054 

; Product MODEM 101 

; Module SYSTEM 

This file -defines the standard systerr definitions. 

Revision Number $Revision$ 

Revision Status- $State$ 

Last Modified: $Date$ 

Identification: $ld$ 



Revision History: SLogS 



7 



#ifndef SYSTEM STANDARD H 
^define" SYSTEM" STANDARtX H 



r 



#ifdefDO CX025 
#inciude <T3LBL Globals h> 
#include <m68Er:040 h> 
#endif /* DO_CX025 7 

#ifdef DO_WIN 

/ 4 -#include < null.h>-7 
#derine STDC FALSE 
//#inclucfe~ <stddeTTi> 

#define strupr(a) __strupr(a) 
#define strlwr(a) _strlwr(a) 
#define inp{aj jnp(a) 
#define inpw(a) mpw(a) 
#define outp(a.bT outp(a.b) 
#define outpw(a.bT _outpw(a.b) 

#ifdefDO WIN32 
#ifdefDO ADI218 1 
#define DEBLEVEL 1 
#define DEBUG 

include <debug h> /*- DDK Debug Header File 

//tfdefine Debug_Printf Service LCODE_Debtjg_Printf_Service 
//#include <vxdwraps h> 
#endif /* 00_ADI2181 7 
#define far 
#define nea r 

#define NULL f(void *)0) 
#endif /* DO WIN32 7 

#else /* DO_WIN 7 

#ifndef DO_CX018 A 
/*#include < null-h^V 
tfifndef NULL 

#if defined(TlNY ) || defined( SMALL_ ) || defined(_MEDIUM_) 

#define NUTT 0 — 

$else 

#define NULL 0L 

#endif 

#endif 

#endif r DO_CX018_A 7 
#endif /* DO.WIN 7 
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typedef unsigned char uin18; /* range 0-255 7 

typedef unsigned char umod8: /* 8 bti math requirements 7 

typedef unsigned char octet: /* octet data 7 

ffirdef DO FORCED BYTE ALIGNED 

typedef uhsiqned shorTpackS; 

#e!se r DO TORCED BYTE_ALIGNED 7 

typedef unsigned char pack8; /* packed octet data 7 
ffendif /* DO FORCED BYTE_ALlGNED 7 
typedef pack"B expand87T expand octet data to uint16 7 
typedef signed char sint7; 

typedef unsigned short uint16. 
typedef signed short sint 1 5, 
typedef signed short sfract15. 

typedef unsigned long uint32. 
typedef S'gned long sint3L 

#ifdef DC_CX025 
#define far 
#define interrupt 
tfendif /* DOJ3X025 7 

/ — — — — * v 

#define MAX UINT8 (255) 
#define MAX "SINT7 (127) 
#define MINJBINT7 ( 128) 

#define MAX UINT16 (65535) 
#define MAX~"S!NT15 (32767) 
#define MINT5INT15 (-32768) 

#define MAX UINT32 (4294967296 
#define MAX~SINT31 (2147483647 
#define MINJSINT31 (-2147483648 

#define MIN DOUBLE r 7e-308) 
#define MAX"_DOUBLE (1 7e308) 



typedef uirM8 byte; 
typedef uint16 word; 

typedef int bool; 

#ifndefDO CX025 
#define FALSE 0 
#define TRUE ('FALSE) 
#endif /* DO_CX025 7 

#if defined(DO ADI2181) || defined(DOJSAR) 

typedef uint32 3sp_pm_t; /* Of size of DSP program memory contents 7 

typedef Hint 16 dsp_dm t; /* Of size of DSP data memory contents 7 

typedef uint16 dsp addr_t; /* Of size of DSP addresses 
typedef uint 16 dsp" xddr t; /* Of size cf extended DSP addresses 7 
#endif /* DO_ADI2l81 || DOJSAR 7 

/ *** * — — 7 

#define forever for (;;) 
#ifdefDO CX025 

#deftne inTerrupt_disable_onto stackf) \ 

asmj "6ri w #0x0700. sr") 
#define interrupt_restore_from stack() \ 

asm( "andi.w #0xF8FF. sr") 

#else r DO_CX025 7 

#ifdefDO CX018 _A 

#define interrupt cfisable onto stackQ \ 

{asrnf ORI #$700,SR\n NORn NOPVn");} 
#define interrupt restore from stack() \ 

{asrnf ANDI #$f8ff.SF«h NOP\n NORn"):} 
#else r DO_CX018_A 7 

#ifdefDO WIN 

#define inferrupt_disable_onto_stack() 
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#define interrupt restore from stack() 
#else/ 4 DOJAJIflTV 

#define interrupt_disable onto_stackQ (asm pushf;asm cli;} 
#define interrupt restorelrom stack() {asm popf;} 
#endif r DO_VV|ft v 

#endif r DO CX0I8 A V 
#endif r DO__CX025" T 7 

#ifdefDO CAPI20 

#define DEBLEVEL 1 

#pragma code_segr_LTEXT". "LCOOE") 

include <detxig.h> 

^include <vxdwraps.h> 

#else r DO CAPI20 V 
#ifdef DO_l3AR 

typedef const char far* LPCSTR; 
typedef char far* LPSTR; 
#define PASCAL pascal 
tfendif /* DO ISAR 7 
#endif /* DO" CAP120 V 

#ifdef DOJA/IN32 
tf#define sprintf Sprintf 
#endif /* DO WIN32 V 

r 7 

r * www... v 

/ — * v 

f * v 

#define print_debug off(A) 
#define print_diag_ oTf(A) 

tfifndef DO_CAPI20 

#define prinl_pr,mary(A) print routine A 

#define print_std(A) print rouFtne A 

#define prinMnfo(A) prinF_rouline A 

#define print_diag_on(A) print rot/tine A 

#define print debug_on(A} print routine A 

#define report anomaly(A) report anomaly routine (A) 

#else /* DO_CAPI20 7 

#define print_std(A) Debug_Printf A 

#define prinMnfo(A) Debug Printf A 

#define print_diag_on(A) Uebug_Printf A 

#define print debug_on(A^ Debug_Printf A 

#define reportanomaly(A) report_anomaly_routine (A) 

#endif /* OO CAPI20 7 

#ifdef DO_PRIN1 SDISABLED 

#undef print pimary 

#undef prin(_"std 

#undef prints info 

#undef print debug 

#undef reporl_anoniaty 

#define print primary{A) 

#define print^std(A) 

#dcfine prinHnfo(A) 

#define print debug(A) 

#define report anomafy(A) 

#endif r DO_ PR I NTS JIM SABLE D 7 

#ifdefDO CX025 
#undef priht_primarv(A) 
#undef print_sid(A) ' 
#undef print_info(A) 
#undef print debug(A) 
#undef report ahomafy(A) 
#endif /* DO_CX025 7 

r **** * — — v 

#endif /* SYSTEM_STANDARD_H 7 
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